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We describe a nonperturbative method for calculating the QCD vacuum and glueball wave functions, based on 
an eigenvalue equation approach to Hamiltonian lattice gauge theory. Therefore, one can obtain more physical 
information than the conventional simulation methods. For simplicity, we take the 2+1 dimensional U(l) model 
as an example. The generalization of this method to 3+1 dimensional QCD is straightforward. 



1. Wave Functions 

To study more preciesely the glueball prop- 
erties, one should compute not only the spec- 
troscopy, but also the wave functions. Once the 
wave functions are available, the matrix elements 
relevant for the glueball production and decays 
(branch ratios) become calculable. In some sence, 
the wave functions can give more physical 
information than the masses themselves. In this 
paper, we present a method for such a purpose. 

We begin with the vacuum wave function in the 
form 



preserves the continuum limit, where the opera- 
tors R and are expanded in order of graphs 
(clusters): 



i? = i?i + i?2 + • • • , F^ = Ff 



F4 



Here i?i (or F^) is the lowest order term in R (or 
F^), and is chosen to be 

i?i-riGi, Gi = i^[i7p(x) + /i.c.]. 



n >-- 



0>, 



h.c 



where R{U) consists of Wilson loops (clusters) 
m, and the state | > is defined hy Ei | >= 0. 
For the glueballs, we consider here only the anti- 
symmetric (under a plane parity transformation) 
lowest lying excited state, and take the wave func- 
tion as 



with a coefficient ri (and ai ) to be determined. 
Higher order clusters can be produced by solv- 
ing the eigenvalue equations for | > and 
I A > order by order, from which wc obtain 



F' 



(1) 



R(U) 



0>, 



where F^{U) contains various Wilson loops with 
the appropriate symmetry. It has been shown 
that we can establish a truncation scheme, which 



* Supported by National Natural Science Foundation of 
China 



Here the repeated index j implies a summation 
over all the clusters up to some n-th order. In 
previous papers we have shown that this 

method is very efficient in obtaining scaling be- 
havior of physical quantities. First results |^ for 
the glueball masses in QCD have been obtained 
from this method. 
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2. Correlation lengths 

In order to obtain the correlation of the states, 
we first investigate the continuum hmit of the 
clusters. Expanding Up (plaquette), Gj, and Gf 
in order of the lattice spacing a, we have 

Then 

Gj ^1- Aja^e^T^ - Bja^e^TiVl + Vl)T + • • ■ , 
Gf = -X^a^e^T^ 



For illustration and simplicity, we consider here 
a 2+1 dimensional U(l) model. It is well known 
fioll that the theory is confining for all non- 
vanishing coupling constant. When a goes to 
zero, the glueball mass Ma is expected to de- 
crease exponentially as 

M a -^exp ^ , 

9 \ 9 J 

where ci and C2 are some constants. Hence, in 
the scaling region 

Ho' = 2 ln{grjAj) = 2 ln{fioMc'^^) + C2/3, 



-Yja^e^T^{Vl+Vl)T+--- (2) 

where T — J-12 is the field strength tensor, Vi 
and I?2 are the covariant derivatives, and Aj and 
Bj or Xj and Yj are constants, corresponding to 
the cluster Gj or Gf respectively. The long wave 
length vacuum wave function is [^|J4[|^ 

According to Eqs. (|l|) and (||), 

/io — rjAjO^e^, 112 = rjBjO^e^, 

from which we obtain the correlation length be- 
tween the field strengths of the vacuum 

= a^TjBjjTjAj. 

Similarly, the long wavelength anti-symmetric 
glueball operator is 

and 

The correlation length in the anti-symmetric 
glueball state is 

u = <^{uY~ii~x, 

At this point, it should be noted that the results 
above are very general in the sense that the spa- 
cial dimension and the gauge group of the theory 
are not specified. 



M2/ = g'"-(.9 ^^j-Sj) -ln{^i2Mc^ ^) -I- C2/3, 
= ln{PrjB,/r,A,) = ln{e,M^c^') + C2P, 
4' = \^n{g^fjXj) = ^ln{iJL^Mic~^^) + C2/3, 

= \ln{9^f,Y,) = ^ln{^i^Aficf') + C2/3, 

W - ln{Pf,Y,lf,X,) = InieAM^c^^) + C2/?. 

In the contimmm limit, /ig, /i2, ^J^Q , M^i Cm and 
M should be constants, which means the curves 

of ^Q, against (3 will be straight line with 

the expected slope in the scaling region. 

In Figs. 1 and 2, we present the results for 
/iQ and ^2 against (} from the 2nd to 5th order. 
In Fig. 3, we also show the result for /x^/ from 
the 2nd order through 4th order against (3. All 
the curves show a nice exponential behavior and 
a clear trend towards convergence. In particular, 
the slope of is well consistent with the spec- 
trum. For /L(2,/Xo /, M^2' ^'v we expect that 
when the order increases, their slopes would ap- 
proach their correct theoretical values. 

In conclusion, the method described above has 
proved to be very useful for calculating the vac- 
uum and glueball wave functions and correlation 
lengths. The calculations in 3-1-1 dimensional 
QCD is in progress and the results will be re- 
ported elsewhere. 
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Figure 1. fiQf as a function of /? 

Figure 2. as a function of /3 
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Figure 3. /ig / as a function of /? 



